Relativistic Dissipative Hydrodynamic Equations at the Second Order for 
Multi-Component Systems with Multiple Conserved Currents 
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We derive the second order hydrodynamic equations for the relativistic system of multi- 
components with multiple conserved currents by generalizing the Israel-Stewart theory and Grad's 
moment method. We find that, in addition to the conventional moment equations, extra moment 
equations associated with conserved currents should be introduced to consistently match the number 
of equations with that of unknowns and to satisfy the Onsager reciprocal relations. Consistent ex- 
pansion of the entropy current leads to constitutive equations which involve the terms not appearing 
in the original Israel-Stewart theory even in the single component limit. We also find several terms 
which exhibit thermal diffusion such as Soret and Dufour effects. We finally compare our results 
with those of other existing formalisms. 

PACS numbers: 25.75.-q, 25.75.Nq, 12.38.Mh, 12.38.Qk 



I. INTRODUCTION 

Hydrodynamics, which is grounded on conservation 
laws under local equilibrium conditions, is widely used 
in general physics. Its relativistic version taking account 
of irreversible processes was initiated by Eckart [l[ many 
years ago. Later, Landau Q applied relativistic hydro- 
dynamics to multi-particle production in hadron-hadron 
collisions in cosmic ray events. Just after the historical 
work, the applicability of relativistic hydrodynamics by 
Landau was examined in terms of quantum field theories 
3. However, the theories by Eckart [Ij and by Landau Q 
share a common problem that dissipative perturbation 
propagates at infinite speed Q, which is obviously in- 
compatible with the concept of causality in the relativis- 
tic theory. The problem is originated from rather phe- 
nomenological constitutive equations for dissipative cur- 
rents: Instantaneous responses to thermodynamic forces, 
which are usually assumed in non-relativistic irreversible 
processes such as Fourier's law and Newton's law, lack 
the relaxation of the dissipative currents. The theory 
is called first order theory when the dissipative currents 
are proportional to thermodynamic forces since the en- 
tropy current in ideal hydrodynamics is corrected by the 
linear terms of dissipative currents. On the other hand, 
second order theory, in which the entropy current has 
quadratic terms of dissipative currents, leads to the re- 
laxation terms for dissipative currents and therefore can 
satisfy the causality. So far, a wide variety of second 
order theories j6l-l2l| have been proposed as relativistic 
theories of irreversible processes. The expression of con- 
stitutive equations, however, varies among theories. 

In this paper, we derive relativistic dissipative hy- 
drodynamic equations at the second order with multi- 
components as well as multiple conserved currents. One 
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finds several approaches to obtain constitutive equations 
for relativistic systems in the literature. First of all, to 
clarify the aim of this paper, we categorize the systems 
into four classes according to the number of components 
and their interaction: (a) single component with binary 
collisions, (b) single component with elastic and inelas- 
tic collisions, (c) multi-components with binary collisions 
and (d) multi-components with elastic and inelastic col- 
lisions. In all the classes above, one of the macroscopic 
equations to be solved is the energy-momentum conserva- 
tion. The typical situation, which can be found in many 
textbooks of kinetic theory or non-equilibrium statistical 
physics, belongs to the first class (a), in which the number 
of particles is conserved during evolution. In addition to 
the energy-momentum conservation, one needs to solve 
the continuity equation of the number of particles. In 
the second class (b), the number of particles is deter- 
mined locally by temperature and is not necessary con- 
served due to inelastic processes during evolution under 
local thermal and chemical equilibrium. Then, one solves 
the energy-momentum conservation only in this case. In 
the third class (c), the number of each component is 
conserved due to binary collisions. Thus, the number 
conservation for each component as well as the energy- 
momentum conservation are solved simultaneously. In 
the fourth class (d), which we will discuss in this paper, 
the number of each component may not be conserved 
due to inelastic collisions or chemical reactions. However, 
there can exist several conserved numbers due to symme- 
try of Lagrangian under some continuous transformation. 
Instead of the number conservation for each component, 
one needs to solve the continuity equations of conserved 
charges together with energy-momentum conservation. 
Note that the number of components does not need to 
coincide with the number of conserved currents. To our 
best knowledge, no systematic investigation is available 
for the class (d) above even though it is the most impor- 
tant situation in ultra-relativistic systems at very high 
temperature in which particle creation and annihilation 
take place frequently. 
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The importance of relativistic hydrodynamics [22[ 
has been increasing after the discovery of the "perfect 
fluid" quark-gluon plasma (QGP) in Au+Au colhsions 
at y/SNN = 200 GeV at Relativistic Heavy Ion Collider 
(RHIC) in Brookhaven National Laboratory (BNL) 
This fact has been quantified by reproducing, within ideal 
hydrodynamic models p3 - l26j , particle spectra as well as 
elliptic flow coefficients V2 (ellipticity of radial flow in mo- 
mentum space d?!) as functions of centrality, transverse 
momentum and pseudorapidity from experimental data 
[28l - l30j . Since ideal hydrodynamic models are approx- 
imation in the sense that all non-equilibrium processes 
arc omitted, our next step should be to include small 
viscosity to capture the correct physics. The reasonable 
agreement of ideal hydrodynamic results with experimen- 
tal data suggests that the system is not so far from equi- 
librium, i.e., viscous hydrodynamic models can be jus- 
tified for the QGP at RHIC energies. In Large Hadron 
Collider (LHC) experiments which have just begun [31 1 
and are planned to eventually reach y/s^N = 5.5 TeV 
in the heavy ion program, viscous hydrodynamic mod- 
els will become even more important in quantifying the 
properties of the hot QCD matter and examining the ap- 
plicability of hydrodynamic models. It should be noted 
here that hydrodynamics and a hydrodynamic model are 
different concepts and thus are to be distinguished; the 
former is a general macroscopic theory that describes 
strongly-coupled relativistic systems, while the latter is 
a specific model based on hydrodynamics that describes 
the phenomena of interest, namely relativistic heavy ion 
collisions. It is essential that viscous hydrodynamics be 
established before constructing any realistic models for 
heavy ion collisions. 

We aim to develop the formalisms of relativistic dissi- 
pative hydrodynamics for multi-component systems with 
multi-conserved currents by determining the distortion of 
distribution functions and then constraining the constitu- 
tive equations for the dissipative currents. The discussion 
for multi-component systems was not recognized well in 
the context of highly relativistic system where particle 
creation and annihilation take place, but it turned out 
to be far from trivial [3^ . Multi-component hydrody- 
namics is important in developing dissipative hydrody- 
namic models for the hot QCD matter, which also is a 
multi-component system. We also consider systems with 
multiple conserved currents because one should be able to 
introduce more than one conserved charge such as baryon 
number and strangeness to the system. 

This paper is organized as follows. In Sec. HIl we show 
how to consistently formulate the second order consti- 
tutive equations for multi-component systems. The sys- 
tems with multiple conserved currents arc considered. In 
Sec. Illli wc discuss the correspondences between our re- 
sults and several different existing second order equa- 
tions. The conclusion will be given in Sec. IIVI The 
Minkowski metric g'^" = diag(-f-, — , — , — ) and the nat- 
ural unit c = h ~ ks ~ 1 are used throughout this 
paper. 



II. DERIVATION OF SECOND-ORDER 
VISCOUS HYDRODYNAMIC EQUATIONS 

We derive macroscopic dissipative hydrodynamic equa- 
tions in multi-component systems with multiple con- 
served currents by extending the Israel-Stewart second 
order theory Q. We consider the third moment of the 
distribution function /* and constrain its derivative from 
the second law of thermodynamics. We also discuss ad- 
ditional moment equations which do not appear in the 
conventional Israel- Stewart formalism to consistently de- 
scribe multi-component systems. 

Wc find several non-trivialities of multi-component sys- 
tems in the course of the formulation. Firstly, the ther- 
modynamic stability conditions, which ensure that the 
system is in maximum entropy state in terms of dissipa- 
tive currents, have to be employed after the constitutive 
equations for dissipative currents are obtained, because 
the number of the constitutive equations and that of the 
dissipative currents would not match if the conditions 
were considered beforehand. Secondly, the second law of 
thermodynamics requires a specific tensor structure for 
moment expansion of the distortion of the distribution 
Sf" in a multi-component system, which justifies the re- 
sult of Ref . [ill . Thirdly, if the system has conserved cur- 
rents wc need to consider new moment equations to con- 
sistently formulate multi-component relativistic dissipa- 
tive hydrodynamics. These equations also allow us to de- 
termine all the dissipative currents in an arbitrary frame, 
which the conventional Israel-Stewart theory would not 
do. 

The existence of multiple conserved currents also 
brings uncertainties to the conventional Grad's 14- 
moment method since it is no longer applicable when 
more than 14 dissipative currents are present. In this 
paper we consider systems with conservations based on 
quantum numbers, such as baryon number, strangeness 
and isospin. In other words, inelastic scattering and 
chemical interactions are present. It should be noted 
again that the number of conserved currents and that of 
particle species in the system are generally different when 
inelastic processes are present. We propose a generalized 
moment method based on Onsager reciprocal relations 
[33! to describe such systems without ambiguity. 



A. Extended Second Order Israel-Stewart Theory 
for Multi-Component Systems with Multiple 
Conserved Currents 

Wc would like to introduce thermodynamic quantities 
in the tensor decompositions of the energy-momentum 
tensor and the conserved currents. In the ideal relativis- 
tic hydrodynamics the energy-momentum tensor and the 
conserved currents are expressed as 



N'^ - 



CqWU 



(1) 
(2) 
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where the index J ( J = 1, • • • , N) denotes different types 
of conserved currents. Here N is the number of con- 
served currents, is the four velocity normahzed as 
m'^m^j = 1. A^'' = (7^'^ — u^u'^ is the projection opera- 
tor, eo = u^T^"^u„, Pq = -iApi.Tp'"' and njo = u^N^^ 
denote the energy density, the hydrostatic pressure and 
the charge number density of the J-th conserved cur- 
rent, respectively. The number of unknowns here is 
5 -|- TV because one has eo(l unknown), Po(l unknown), 
^jo unknowns) and u^(3 imknowns). On the other 
hand there are only 4 -I- iV equations from the energy- 
momentum conservation and the charge number conser- 
vations to describe the system: 

9,r^'' ^ d^T^" = 0, (3) 
d.N^j = d,N% = 0. (4) 

Therefore, one needs to introduce the equation of state 
= ^o(eo, {'T-jo}) from microscopic physics to com- 
pletely determine the space-time evolution of the system. 

On the other hand, the energy-momentum tensor and 
the conserved currents in relativistic dissipative hydro- 
dynamics are tensor-decomposed into 

T^"^ = {ea + 5e)u^'u'' - {Po + I{)l\^"' 

+ 2W^(^w'') +7r^'', (5) 
N'^j = {nj^ + 5n.,)u^^ + Vlf. (6) 

We denote the dissipative parts of the above equations 
as 5T^''' = T'"' - T^" and 6N^ = - ^,/o- Then 
n = -^A^JT'"' is the bulk pressure, = A^'JT°'^up 
the energy current, tt'"^ = (5r<^'^> = [i(A^„A''^ + 
A^^A''^) - lA^"'Aa,p]ST°'^ the shear stress tensor and 
Vj = Af^^SN"^ the charge current of the J-th con- 
served current. Round bracket for Lorcntz indices de- 
notes symmetrization as A^t^B"^ = ^{A'^B" + B^'A"). 
5e = Uf^ST^^^Ui, and Snj = Uf^SNj are the distortion of 
the energy density and of the J-th charge density, respec- 
tively. They actually vanish because the stability condi- 
tions need to be employed. The details can be found in 
Appendix of Ref. [321 . However, we have to keep these 
quantities for the moment to correctly count the number 
of unknowns. 

The energy-momentum conservation and the charge 
number conservations in non-equilibrium systems are 
tensor-decomposed into 

D{eo + 6e) = -{eo + 6e + Pa ^u^' 

+ 2W''Du^-V^W'' 

+ ^^"W^^u,^, (7) 

{ea+5e + Po + \l)Du^' = V [P^ + \1) - ^u" 

- Af^DW^-WW^ui" 

+ ^^''I^ii,, - A^-'V^.p, (8) 

D{njo + Snj) = -(n,/o + 5nj)V ^u^ 

' V^V^ + V;;Du^, (9) 

where the time-like and the space-like derivatives are de- 
fined as D = and = A^^d" , respectively. 



There are 10-1- AN additional macroscopic variables 
H, 5e, W^, TT^'', 6nj and V^, which originate from 10 
independent components of ST'^'^ and AN independent 
components of SNj. Therefore one needs to introduce 
constitutive equations for these dissipative currents to 
completely describe the system. The relativistic Navier- 
Stokes equations of the dissipative currents can be ob- 
tained by considering the law of increasing entropy. The 
energy momentum tensor T'"^, the conserved currents 
Nj and the entropy current arc expressed in kinetic 
theory with the microscopic phase-space distribution 
as 

where gi is the degeneracy and the conserved charge 
number of the J-th conserved current. 0(/') = /* In /* — 
e~^(l -I- e/*) ln(l + e/*) where the sign factor e is +1 for 
bosons, —1 for fermions and for classical particles. In 
the classical limit, </)(/*) reduces to /* In/* — /*. 

We write the expansion of the entropy current up to 
the second order in 5f ^ ~ /' — /q as 

s'^ = s'^ + 5sl^^+5s^l^^+0{5f). (13) 

If one naively considers the first order dissipative cor- 
rection then one obtains the Navier-Stokes expressions 
of the relativistic constitutive equations from the law 
of increasing entropy. Here we would like to empha- 
size that as is the case for standard statistical mechanics, 
one needs to have Onsager cross terms in the first order 
expressions, even though they are sometimes neglected 
in conventional relativistic formalisms. The cross terms 
are actually essential because they give rise to impor- 
tant physics such as Soret and Dufour effects and also 
play a significant role in preserving the law of increas- 
ing entropy. Detailed discussion on the relativistic linear 
response theory can be found in Appendix [X] However, 
the formalism is known to be both acausal and unstable 
Q as they allow propagation of information faster than 
the speed of light. To avoid such unnatural behavior, one 
has to consider the second order correction to s^ which 
introduces relaxation effects on the dissipative currents. 

Israel-Stewart formalism of the second order single- 
component dissipative hydrodynamics can be derived 
from the third moment of the distribution /'. We gen- 
eralize the formalism to multi-component systems and 
write 

a./'^-" = ^ I -^^ptp'ip-d^r = (14) 

i 

Here F^'' is a symmetric tensor which we will deter- 
mine later. There are 10 independent equations, which 
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are, when projected either parallel/perpendicular to the 



flow u'^, the scalar equations u^j^Uudal'^'^ 



(1 equation) and /^^^d^I'"'"' = A^^.F^'^ (1 equation), 
the vector equation Apfj^u^dal^"'^ = Ap^u,yF^'' (3 equa- 
tions) and the tensor equation 9^/^'"'^" = y(A'i'> (5 equa- 
tions). 

On the other hand, the number of unknown dissipativc 
currents is 10 -|- AN in the systems with non- vanishing 
chemical potentials. In this respect, we need AN more 
equations to fully determine the system. Since the con- 
ventional moment equations are the derivatives of the 
higher order moment for the energy-momentum conser- 
vation as shown in Eq. (jl4p . it would be natural to in- 
troduce new moment equations in multi-component sys- 
tems, 



Ualj — 



E 



(15) 



which are the derivatives of the higher order moment 
for the charge number conservations. One might argue 
that weight factors other than the conserved charge num- 
ber qf could have been chosen to construct the moment 
equations, but Eq. (jlSp is the only vector equation which 
vanishes in the limit of zero net charge number densi- 
ties. This formalism not only makes the correspondence 
between the number of equations and that of unknowns 
clear, but also allows one to construct an arbitrary num- 
ber of moment equations of the form dalj" = Yj de- 
pending on how many conserved charge currents the sys- 
tem possesses. 

We now have 10-|-4A^ equations for 10-|-4A^ dissipative 
currents. The next task is to estimate and Yj to de- 
rive the constitutive equations in terms of the dissipative 
currents. The constraints arc again given by the second 
law of thermodynamics. The definition of the entropy 
current p2|) gives 



-E 
E 



ptv'd.r 



where is defined in = [exp(y') — e]^^. As discussed 
in Appendix the second order constitutive equations 
depend on an explicit form of the off-equilibrium distri- 
bution /% or equivalently, y'. When the deviation from 
the local thermal equilibrium is small, we may write 
as 

f^vl + 5y\ (17) 
where is defined in /q = [cxp (t/g) — e]^^, which means 



E 



(18) 



fij is the chemical potential associated with the J-th con- 
served quantity. We estimate the off-equilibrium correc- 
tion dy^ through Grad's moment method. The conven- 
tional Grad's 14-moment method 0,113] cannot be gener- 
ally applied to the systems with N conserved charge cur- 
rents because the number of macroscopic variables is then 
10-1- AN whereas that of unknowns in the expansion re- 
mains 14. Actually, if the system has no conserved charge 
currents, i.e., N = 0, the system also cannot be solved, 
unless the concept of one conserved current is present in 
the system. This means that the system of single com- 
ponent with binary collisions or multi-components with 
one conserved charge current is implicitly assumed in the 
Grad's 14-moment method. In the latter case, if there is 
no conserved charge current, the limit of vanishing chem- 
ical potential for the stability condition can be taken js^ . 
Thus we have to introduce 10 -I- 4A^ unknowns in the ex- 
pansion as well. If we demand that (i) the distortion of 
the distribution can be expressed in terms of the dissi- 
pativc currents, (ii) the resulting constitutive equations 
for the dissipative currents satisfy the Onsager reciprocal 
relations and (iii) an arbitrary number of conserved cur- 
rents can be introduced to the system and the effects of 
a conserved current vanish in the vanishing limit of the 
corresponding chemical potential, then the only possible 
way is to assume the expansion 



6f = 



(19) 



where and e^^ are macroscopic coefficients of the ex- 
pansion which include information of all the components 
in the system. We will see later that this indeed yields the 
second order constitutive equations which are reasonable 
compared with those of other formalisms. 

We make several comments on this distortion Sp. 
Firstly, a non-zero trace tensor correction e^^i, is con- 
sidered instead of a scalar and traceless tensor correc- 
tion in Eq. p^ . We will see later that this must be 



the case for multi-component systems [32[ because of 
the law of increasing entropy. Secondly, if one consid- 
ers a single component system with binary collisions, 
i.e., q = 1, then Eq. ([T9| reduces to the expansion 
(16) Sy = p^Efj, +p^p^e^i, which is equivalent to the conven- 
tional Grad's 14-moment method mentioned in Ref. 0]. 
Note that this is not completely equivalent to the Grad's 
14-moment method for the systems with chemical inter- 
action. Thirdly, the distribution function satisfies the 
stability conditions because 6e and 5nj are treated care- 
fully in this formalism and they can be considered zero 
after the constitutive equations are derived. 

The unknowns in the distortion of the distribution can 
be determined by matching the macroscopic variables 
and the ones calculated in relativistic kinetic theory. We 



first tensor-decompose the 10-f 4A^ unknowns e^i, and 
in terms of the flow as 



A,, 



n Ml 



(20) 
(21) 
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where we employed the notcitions — ^fiv^ £\\ 



J 



expressed as 



A'^'ei and 
Then the consistency conditions can be 



Jo 



V ••• J 



( 



Ji 



V ••• J 

J2£ I I 



V ••• J 

V ••• / 



where the matching matrices are defined as 



/ |J42 J41 Ji\ J3I 



Jo 



J41 Jdi 



^31 



'3 
'^30 



tJ2 
'^30 



'^30 '^20 '^20 



p2 T. 

■^31 '^30 



32 J]l]2 J3232 

•^20 "^20 



Jl 



J2 



32 

( 2^/41 Ji\ Ji\ 

9 tJI jhil jjlj2 

^"^31 "^21 ''21 

'32 

31 



V 

-2J4: 



9 t32 t3i32 J3232 

^Jcii 0/21 ■J2I 



J 



Here the moments of the distribution are defined as 



E 



n 

+ (permutations)] Jjj. 



/o(i + 



,M2„ + 1 . . . , 



jfc- 



(25) 

(26) 
(27) 

) 

(28) 



where the index jk ■ ■ ■ denotes the additional weight fac- 
tor Qj'qf ■ ■ ■ in the summation over particle species i. 
The moments with no index mean no weight factor in the 
summation. Note that the moments with double charge 
weight, e.g., q^q^ do not vanish in the limit of vanishing 
chemical potential. 

Then Eqs. can be easily solved and we obtain 



J 

(Bull + BseSe + ^ I3snjSnj)u^u„ 
J 

2BwU(^W,) +2j2BvjUipVj^ + B^^^„(29) 



,7 



K 



K 



(30) 



where _B's and ZJ's are macroscopic quantities which can 
be determined by temperature and chemical potential 
only and are momentum independent. They contain in- 
formation of all the constituent particles in the system. 
The explicit expressions of the prefactors for the dissipa- 
tive currents in relativistic kinetic theory are 



(22) 


Bn = 


~g(^0 Bse — g(Jo ^)l2i 


(31) 




Bsnj = 


g(Jo"^)i,2+j, -Bn = -(Jcr^)2i, 


(32) 


(23) 


Bse = 


(Jo ^)22, Bsnj = (Jo ^)2,2+i, 


(33) 


Bw = 


(Jf ^)ll, Bvj = (Jf ^)i,i+j, 


(34) 




BjT ~ 


(J2 )ll, -Djl = -(Jo )2+j,l, 


(35) 


(24) 




(Jo ^)2+j\2, Dj„^ = (Jo ^)2+J,2+fc, 


(36) 






(Jl ^)l+i4' D'yk = ^Jl 


(37) 



where 1 < (fc, /) < N . It should be noted that 5y^ is ex- 
pressed as linear combinations of the dissipative currents. 
Higher order contributions are discussed in Appendix |B] 

The entropy production p6| is now expressed as 



E 

■i 

E 



3p.(2/o + '^2/^K5^r 



{2^fE, 



E^) 



(38) 



where we have used the energy-momentum conservation 
and the charge number conservations. The tensor struc- 
tures of £vp and ell in terms of the flow needs to 
be reflected on Y^f and YJ in linear response theory. 
Note that again the "cross terms" are allowed for the 
reason mentioned in Appendix [X] The finite trace tensor 
correction to the distribution e^^ is required instead of 
the tracelcss one with the scalar correction e for systems 
with multi-components and/or multiple conserved cur- 
rents, because 7^ can be shown from the fact that 
the trace of Eq. ([H]) is not generally zero, i.e.. 



{2TTfE, 



■py^r ^ 0, (39) 



K 



and the existence of the moment equation for e violates 
the matching of the number of dissipative currents and 
that of the constitutive equations. It gives a justification 
to the expansion Eq. ([T^ that a non-zero trace e'^" should 
be considered. The tensor structures of Y'^'^ and Yj are 
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then expressed as 

( ^r^ \ 



= Cn 



(40) 



Y 



V ••• / 



V ••• / 



where 



Y,,,. 



Y\\\ 



Ui_,u^ + Y 



2Y,"-u 



Y. 



Y 



(41) 
(42) 

(43) 
(44) 



The microscopic physics is integrated out in the trans- 
port coefhcient matrices C^'s. Here C^'s are semi-positive 
definite and symmetric because of Onsager reciprocal re- 
lations. Yj — and — )• in the limit of corresponding 
vanishing chemical potential /ij implies that the trans- 
port coefficients for the cross terms between and Y^l ^ 
or equivalently, and Y^^i,^ also vanish in the limit. We 
uniquely determine the constitutive equations for all the 
dissipative currents by solving Eqs. p0|) - (|42l) in terms of 
the 10 + 4iV dissipative currents e^^ and ej and then in 
terms of 11, 8e. W^, tt'^", Snj and Vj using the explicit 
forms of the distortion of the phase-space distribution 

We obtain the constitutive equations explicitly by es- 
timating the derivatives dal^'^°' and dal'j" in Eqs. (fT4|) 

Remembering y' = ^^T^ + Pf("T" + 



and lfT5|). 

order are 



their expressions up to the second 



T 



UaCR 



E 

J,K 



J 



K 

K 



da 



T 



T 



(45) 



8 T'^" 



K 

E^^^^^'+E^r'^O^-f 

- r T / 



where the additional moments are defined as 



^jk- 



E 



X (l + 2e/^)pf^ 



■n 



(permutations)] 



k— 

ran 



(47) 



Again the index jfc • • • denotes the additional weight fac- 
tors g/f?!^ • • • in the summation over particle species. The 
terms in the third and fourth lines of Eq. (j45|) which 
involve the moments if s arc assumed to be small and 
simply omitted in the Israel-Stewart formalism Q. It 
is argued in that paper that the terms proportional to 
the Navier-Stokes thermodynamic forces, i.e., D^, D^, 
V^u'', V^^, and V^^Mi,) may not be significant, 

even though the terms proportional to acceleration Du'^ 
are kept in their equations. However, these are actu- 
ally of the same order in dissipative currents as the other 
second-order terms and should not be neglected to pre- 
serve consistency. Equation ()46p corresponds to the new 
moment equations which do not appear in the original 
Israel-Stewart theory. 

It is now a straight-forward task to derive the constitu- 
tive equations in multi-component systems with multiple 
conserved currents. We have from Eqs. (fT4|) . (fT5|) . ([22])- 

m and mh-imi, 



Sni 
Sn2 

V ••• J 



dal. 
del 



(48) 



/ 2a„/-Lii^" \ 



= -JiCi 



V 



'J a ^ J, 



J2C2daI 



(49) 
(50) 



K 



e,5]da^, (46) 



and when combined with Eqs. (j45|) and (j46|) . the second- 
order constitutive equations for the dissipative currents 
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are expressed as follows: 



n = -CV^u^-rni^n 
J 

J 

J,K r 



(51) 



T 

K.L 



K K 



K 



K 



K 
b 



+ xi.wWV.u^^ + Xv.wWV^u, 

K 

+ Exv^'nnv^^ + xUnv^ 



AT 



(53) 



J 



J 



J J 

E^w.^;^^ 

ExW.^1'^^ +E^W,^7"v''"- 

.7 ,/ 



E^w,,^;v.^^+ExW; 



Ex^i^.TT^'^V.^+xty.'r'^'^V. 



T 



Ex^'nnv^^ + xWnv'l 

,7 

xWn^w^ + xwnV^n, 



1 



Ex^ivT^'^^^+xL^^'-i^;^ 



+ x^.7r^'^Vpu'' + x$J.7r''<^Vpii^ 

1 



+ xUW^^^^^''^+xVV<^VF''> 

+ Ex$J'(:.^i'^v^^^ + Ex^K,^i''V''^^ 

7,74" J 



Ex^v.W"^"^^+E^'^.^^"^;^ 



X.nnV<^7/'^>. 



(54) 



Here C's, k's and r] arc the first order transport coeffi- 
cients which are expressed in terms of the matching ma- 
trices Ji and the semi-positive definite matrices C^. As 
we will see in Sec. Ill Bl these transport coefficients satisfy 
Onsager reciprocal relations. C is called bulk viscosity, 
energy conductivity, r\ shear viscosity and kv> charge 
conductivity of the J-th conserved current, r's arc the 
relaxation times and x's are the second order transport 
coefficients. The stability conditions arc employed at this 
point to take out the constitutive equations for be and 
(577j, and to obtain the second order constitutive equa- 
(52) tions for 11, W^, Vj and ir^" . Apparently we also have 
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the term x^^Tr^'^^V^-Up in Eq. (O, but this term actu- 
ally vanishes because orthogonality relation tt^ 



= 

demands x^tt = 0- We have utilized the fact that the 
prefactors D's and B's appearing in d^Sf^ are functions 
of J^'^ 's, and that their derivatives can be expressed in 
terms of the Navier-Stokes thermodynamic forces: 



K,L,---m,n ""^mn 

V T V FJ^'-d ^ 

"S^ \^ R-fc'- - P) i 



X Jkl - 

K.L.---m,n ^''nin 



(55) 



Here F denotes the prefactors B's and D's. Note that if 
we define the symmetric traceless thermodynamic force 



Vw^). then TT^'^'^V, 



and the vorticity 
in Eq. is 



expressed as 



which is also a commonly found expression. Likewise, 
the identities 

W^l^u^ = ly^'a''^ + VK'^cj^^ + ivP^^V'M^, (57) 
yjV^^.. = y/a^, + VOW', + iv^l-V^w., (58) 

are often used. One should be careful that uj^^ is an anti- 
symmetric tensor. In the first order limit, these consti- 
tutive equations reduce to the Navier-Stokes forms. 

Our formalism have four major differences from the 
conventional method by Israel and Stewart Firstly, 
several second order terms which do not appear in the 
Israel- Stewart theory can be found in Eqs. ([5T|) - ([M1) . 
Actually, these terms - the terms composed of a dissi- 
pative current and a Navier-Stokes thermodynamic force 
- are the results of consistent expansion and should also 
exist in the case of single component systems. We will 
compare our formalism with others in Sec. IIIII and see 
most of the second order terms reported in other papers 
are found in our formalism. Secondly, we can now calcu- 
late the independent second order equations for the vec- 
tor dissipative currents and V^, which allows us to 
determine these variables in an arbitrary frame, whereas 
the Israel- Stewart method yields the three equations for 
q'' only. Thirdly, we have different kinetic expressions 
of the transport coefficients C's, k's, yy, r's and x's due 
to the new moment equations. Note that once the first 
order transport coefficients arc given, one can estimate 
the second order ones and the relaxation times because 
they are related within the framework of kinetic theory. 
Fourthly, chemically interacting systems with multiple 
conserved currents can now be uniquely determined. 



B. Onsager Reciprocal Relations 

We investigate the extended Israel- Stewart theory in 
Sec. Ill Al keeping terms up to the first order and see 
that Onsager reciprocal relations are indeed satisfied in 
our formalism. The entropy production is, according to 
Eq. p6)) . expressed up to the first order as 



p';Sy^d,fl, + 0[d{6f)] 



(59) 



Then semi-positive dcfinitcncss of the above equation 
yields 10 -|- AN constitutive equations. It is straight for- 
ward to derive the first order constitutive equations as 



' 6e * 

Srii 

V ••• J 



-Df 



\ 




\ 



- A 



T 



(60) 



(61) 



(62) 



The explicit forms of the transport coefficient matrices 
are expressed as 



Aq 
Al 

A2 



using the moment equations ([14 
dissipative currents with e's 
thermodynamics and 



(63) 
(64) 
(65) 

T5|) . the matching of 
, the second law of 



illl" 



V 



dalo 

d 



— Jo 



-Df 



(66) 



M T 
-V ^ 



J2 



T 



(67) 



(68) 



Since Ci, the transport coefficient matrices for e^'^ and 
Ej, are symmetric, Ai^i are also completely symmet- 
ric and Onsager reciprocal relations are satisfied. Also, 
^'s are semi-positive definite because C's are. These are 
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the linear transformations of the dissipative currents and 
of the thermodynamic forces in the entropy production 
mentioned in Appendix \^ Note that if the distortion 
of the distribution other than Eq. were employed, 
the reciprocal relations would not hold since the moment 
equations would no longer be uniquely constrained from 
the second law of thermodynamics. Here we emphasize 
that real hydrodynamic transport coefficients reflect mi- 
croscopic physics of the dense medium and are different 
from the ones obtained in kinetic theory and that cal- 
culation of transport coefficients is not the aim of the 
present paper. 

It should be noted here that the derivatives of mo- 
ments in Eq. (j59p . Sq/q"" and daljQ, do not disappear 
because /o is the distribution for local thermal equilib- 
rium, not global one. These derivatives, as we have seen, 
are the sources of Navier-Stokes thermodynamic forces. 
The dissipative currents disappear in local thermal equi- 
librium because the transport coefficients vanish, not the 
thermodynamic forces. 



,l^K 



T 



K=tJ 



K.L 



K 



K 



K 



+ Exv>''''v.^ + x^,.-^''v4 

K 

+ E<nnv^^ + xt,,nnv^i 



K 



T 



(70) 



C. Energy and Particle Frames 



The constitutive equations for multi-component sys- 
tems with multiple conserved currents (j5ip - (|54p are 
frame independent. On the other hand, practically 
speaking, it is convenient to simplify the constitutive 
equations without losing generality by choosing frames. 
There are two conventional ways of choosing a frame 
in dissipative hydrodynamics: the energy frame and the 
particle frame. They are also known as the Landau frame 
and the Eckart frame, respectively. In the energy frame. 



we set the flow u'^ 



in the direction of the flow of 



energy so that no leak from fluid elements exists, i.e., 
Wf' = 0. Note that this implies T'^^w^ = cqu'^. The 
constitutive equations then reduce to 



+ E Xnk^D'-^ + Xnn^D^ + XnnHV^^ 



CnSeD— + E CnSnjD 



T 



J,K 



+ XHttTT^i^V^^U^^ 



(69) 



J 

+ x^.7r^'^Vp< + x:J.^''<^V,4 

+ Ex^'i-.^i'^^^^+Ex^v^.^i'^'^^T 



.J,K 



T 



J J 
+ x.nnV<^4\ (71) 

Note here that the term kvjw{^Du% + \7''^) in Eq. ((70)) 
does not vanish even though it contains the first or- 
der thermodynamic force for W^, because the consti- 
tutive equations for include the term proportional 
to KvKVjV''^ in turn. These phenomena are known 
as Soret effect and Dufour effect respectively in non- 
equilibrium statistical mechanics [ssj . 

The particle frame in single conserved current systems 
is defined as the frame where no leak of the charge is ob- 
served. Naively this is not well defined when more than 
one conserved current is present because in the frame 
where Vj = 0, the other currents would not vanish, i.e., 
^ for J 7^ K. In this respect we should consider 
the average particle frame where the sum of the charge 
dissipation vanishes in the case of multi-conserved cur- 
rent systems. We define here the heat current in a 
system with N conserved currents as 



E 



Pn 



(72) 



which corresponds to the energy conduction through pure 
heat conduction because the contributions from particle 
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diffusions are subtracted. In the average particle frame 

one expects qf^ = W^, or equivalent ly, ^ = 0. The 
flow in this frame can be written as 



''N 



(73) 



which is the average of the flows in one-current particle 
frames. The resulting constitutive equations have the 
same tensor structure as shown in Eqs. ([5T|) -([M 1) because 
each charge current Vj does not vanish. 



III. DISCUSSION 

Several comments arc in order here. 

Firstly, we have derived the constitutive equations 
from the law of increasing entropy only. Actually this 
should be the case for dissipative hydrodynamic formal- 
ism because it is the only thermodynamic relation which 
implies irreversible processes. This is in good contrast 
to the fact that ideal hydrodynamic equations of mo- 
tion, which describe time reversible processes, are the 
energy-momentum conservation and the charge number 
conservations. 

Secondly, if one decomposed ST^"^ and SNj into each 
component i, one would have obtained the equations 
which might be much similar to the single component 
constitutive equations. The problems arc, however, that 
(a) some transport coefficients such as bulk viscosity can- 
not be trivially separated into components, (b) in ad- 
dition to the constitutive equations one needs to solve 
the energy-momentum conservation and the charge num- 
ber conservations but they do not hold for each particle 
species and (c) the distortion of the distribution cannot 
be determined in such methods without introducing addi- 
tional microscopic physics, which often causes lack of gen- 
erality. Such constitutive equations would not be equiva- 
lent to the constitutive equations we obtained. Moreover, 
such equations can be applied only for the system with 
(quasi-)particle picture. In general, it is not the case, in 
particular, in the vicinity of phase transition or in highly 
dense system. 

In the following we would like to investigate other ap- 
proaches and discuss the correspondences between these 
approaches and our formalism. 



requires the information of Sf^. If we alternatively ex- 
pand Eq. (|A2[) up to the second order, we obtain 



E 



-/^(l±/^)(l±2/^)5y'Xyo 



(74) 



which is equivalent to Eq. ()16|) at this order. In this 
case the last term is problematic in obtaining second 
order equations; it involves the terms with two dissipa- 
tive currents coupled with one thermodynamic force, e.g., 
nVF'^V^y when 6y^ is estimated in the moment expan- 
sion. These terms cannot be naively associated with one 
of the dissipative currents to forcefully obtain 10 + 4iV 
equations because generally the dissipative currents of 
different tensor structure can be found in the second or- 
der terms, e.g., V^II terms in the equation for and 

We further consider whether a more phcnomcnologi- 
cal approach 0, [31 ^ in which one expands the entropy 
current with respect to dissipative currents and uses the 
second law of thermodynamics, yields full second order 
constitutive equations. If one assumes that the second 
order distortion of entropy current (jA12[) can be naively 
written as the sum of all the possible second order terms 
in the dissipative currents, one has 



J.K 

+ a^^'USe + a™" 'n<5nj + ^ a^"^'" ' Se6nj 
J J 

J 

+ E ^o'^^'^VJVj' + aJV^'^^.p)"'^ 

J.K 

J 

+ J2 c^i"""-' SeVj^ + J2 ai'^'^SnjW^ 
J J 

J.K 

+ Y.aX-'^Vji:>^\ (75) 



A. Ambiguities of Second Order Equations in 
Phenomenological Approaches 

We calculate the entropy production up to the sec- 
ond order and investigate the possibility of deriving 
10 -|- 4A^ second order equations from the law of increas- 
ing entropy by extending the approach mentioned in Ap- 
pendix]^ We find that second-order equations cannot be 
uniquely determined in this way. As mentioned before, 
the derivation of the second order constitutive equations 



where a's are coefRcients. Again the derivative of the 
entropy current involves the terms with two dissipative 
currents coupled with one thermodynamic gradient, so 
we cannot determine the second order constitutive equa- 
tions in this way. Such terms are naively dropped in 
Ref. 0. 

It is note-worthy that if one considers kinetic theory 
with our extended Grad's moment method to estimate 
Eq. (|A12p . the resulting entropy current will have the 
same tensor structure as Eq. ((75|) does. In this case 
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the coefficients a's are fixed in kinetic theory. When 
one writes down the entropy current first then takes its 
derivative this way, however, the terms proportional to 
Htt'''^ do not seem to appear in s'^ and consequently in the 
constitutive equations since the only possible way to con- 
struct Lorentz vector from bulk-shear term at the second 
order is lin^^Ui, = 0. This does not contradict our results 
because there remains the ambiguity that we can add ar- 
bitrary amount of n7r^''9^u,y — 7r^'TIc)^Uy(= 0) to the en- 
tropy production, and associate one term with the consti- 
tutive equation for the bulk pressure and the other with 
the ones for the shear stress tensor. This corresponds to 
formally keeping a^'^Iin^^'^ u^(= 0) in the entropy cur- 
rent. In other words, constitutive equations from kinetic 
theory in general can have the bulk-shear terms, although 
most of the conventional formalisms seem to be unaware 
of it, possibly because they cannot determine the amount 
of such terms. 

Thus the constitutive equations cannot be uniquely 
determined by the phenomenological approach, nor by 
simply taking derivative of the expansion of the entropy 
current (|A12[) in kinetic theory. The ambiguity of as- 
sociating terms to constitutive equations, of course, is 
removed in our extended Israel-Stewart formalism. 

B. Single-Component Systems without Chemical 
Interaction 

In the conventional second order theories, includ- 
ing Israel-Stewart theory, only the moment equations 
dal'^'^" = Y^'^ are considered. Naively speaking 

the number of equations is 10 in this case, but if the 
single component systems with no particle production 
nor annihilation are assumed the particle number cur- 
rent iV becomes the conserved current, and the num- 
ber of equations is reduced to 9. This is because the 
trace of the moment equations .g^i/^a/^"^" = giiuY^" co- 
incides with the number conservation dfj,N^ = 0, i.e., 
i9q/^^" = rri^dfj^N^ = in kinetic theory. The two scalar 
equations Uf.u^dal^"'"' Uf,u^Y^"' and A^^dal^""" 
Ap^y^'' then become identical. On the other hand, the 
number of dissipative currents is 14 but can also be re- 
duced to 9 by using the stability conditions Se = 5n = 
and choosing the frame to drop 3 dissipative currents 
from either or V^. This means 11, and tt^'' are 
considered. Here the heat current in single component 
systems is defined as 

qt- = Wt'- + (76) 
no 

which reduces to either or y depending on whether 
one chooses the particle frame or the energy frame. How- 
ever, the number of equations and that of dissipative 
currents no longer match in multi-component systems, 
because as mentioned earlier 

g.ud^i'""' ^Y."^lJ j^^^pf^o^r ^ 0, (77) 



is not a conserved current even if no particle production 
nor annihilation is assumed. In a system with particle 
creations and annihilations which we consider in this pa- 
per, what one really has is the charge number conserva- 
tions d^Nj — 0. Thus we have 10 moment equations in 
multi-component systems. This means that naive gen- 
eralization of the conventional second order theory to 
multi-component systems does not work because only 9 
dissipative currents are considered there. Also if the sys- 
tem has more than one conserved current, the number of 
vector dissipative currents exceeds that of equations even 
if the frame is fixed. 

The apparent inconsistency between the number of 
equations and that of unknowns arise from the three 
facts. Firstly, the stability conditions Se = Snj = are 
employed prematurely and 1 -|- unknowns are omitted. 
The stability conditions are employed to ensure that the 
system is in maximum entropy state, i.e., thermodynami- 
cally stable, and thus are different physics from the count 
of the number of unknowns in kinetic theory. In this 
sense the stability conditions have to be considered after 
the constitutive equations are derived. Secondly, the mo- 
ment equations 9q/j" = Yj ([T5|) are not considered for 
the systems with conserved currents. This means that 
AN equations are missing in the formalism. Thirdly, in 
the conventional approach, one considers the heat current 

instead of the energy current and the charge cur- 
rents Vj , possibly due to the limitation of the number of 
equations. This means 3 out of 3 -I- vector dissipative 
currents are taken into account. Thus it is impossible to 
determine all the vector dissipative currents and Vj 
simultaneously in an arbitrary frame within the conven- 
tional Israel-Stewart approach. Related discussion can 
be found in Appendix [Cl 



C. Correspondences with Other Formalisms 

We would like to discuss correspondences between our 
formalism and other frameworks. It should be noted that 
naive comparisons cannot be made because all the other 
formalisms are derived in systems with single component 
and/or single conserved current and are sometimes frame 
dependent. Therefore we take specific conditions such as 
the energy frame or the particle frame in the single con- 
served current limit to make the correspondences clear. 



1. Constitutive Equations from AdS/CFT 

We compare our results with the conformal equations 
for TT^'' based on Anti de-Sitter Space/Conformal Field 
Theory (AdS/CFT) correspondence [1^. Our constitu- 
tive equations for the shear stress tensor can be ex- 
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pressed as, in the conformal limit 11 = 0, 



~ Atttt" 



(78) 



when estimated in the energy frame and in the zero net 
charge density Umits. The identity ([55|) is used here. The 
first term is the Navier-Stokes term, the following two 
terms are the conventional second order terms, and the 
terms in the second line are the new terms that appear 
when the consistent expansion is performed as mentioned 
in the previous section. It is note-worthy that the term 
proportional to IIV omitted in the conformal limit 

but is equivalent to the n^^^V pU^ term in second order 
theory, because both reduces to — C(or 27;)V 
when the first order expressions are utilized. The former 
vanishes while the latter does not, because the transport 
coefficient Q vanishes in the conformal limit, not the ther- 
modynamic force itself. 

On the other hand, the constitutive equations from 
Ref. [ll] in flat space are 

+ ^ttP^Ptt^I - ^TrP^t^Lu'^l + X3LjP''Plj''1 (79) 

where notations are adjusted to our formalism. This 
holds for single component systems, but direct compar- 
ison with our formalism can be made because the con- 
stitutive equations of the shear stress tensor is free from 
the non-trivialities of multi-component systems. Note 
that the terms proportional to tt^-'^D^ in Eq. ([75| can 
be absorbed in the ttP-'^V puP term, because of the ideal 
hydrodynamic relation 



TD- 



1 

T 



(80) 



Then obvious correspondences can be found for all the 
terms, except for the term; the terms with no 

dissipative currents generally do not appear in our for- 
malism, because the distribution is expanded in terms 
of the dissipative currents with Grad's moment method. 
The AdS/CFT approach yields that term because 
all the possible terms which are consistent with their 
approach are added manually in the derivation. If one 
added the vorticity-vorticity term phenomenologically in 
our formalism, that would be inconsistent in the context 
of the derivation in kinetic theory. 



2. Constitutive Equations from Renormalization Group 
Method 

Next we investigate the correspondences between the 
constitutive equations obtained in this paper and the 



ones from renormalization group approach |la | . Here we 
consider a single conserved current system. In the energy 
frame the constitutive equations from that paper are 

n = -CV^uP-mDU 



+ m 
1 

+ 2 
+ Inv 
— Inv 



1 TC „ /■ mu 
2 Tn 



TC 



D 



T 



T 



n 



T5 



nv 



1 



—Du,, 

rp M 



VP 



(81) 



yp 



^) VP^-ryAP^DV^ 



TV 



1 K 



2 TV \ eo + Pq 

r(0) 



9. 



Tyu" f eo + Pq 



noT 



l[-D^+TSi"^D^ + ls^^^V^un^VP 



+ Ivtt 



-V,^+TSv.(v,^ + ^Du, 



+ I 



vn 



-"^P^+TS 



T 



vn 



VP- + -DuP 
T T 



n 



(82) 



ttP" = 2i^V^Pu'''> - t^Dtt^'P-''' 



+ Tt, 



1 Tf] 



2 \ Ttj 
T 

+ T^Si^Ui^P^Pa^l 



- Uv^^PV""! +Ui{V^Pu''^Vl 



(83) 



Our formalism includes the vorticity terms loP'^Vu and 
T^P^Pj^l whereas Eqs. and ([Ml) do not. On the other 
hand the terms which involve derivatives of transport 
coefficients do not exist in our formalism, because such 
terms cannot be expressed with dissipative currents, un- 
less we assume kinetic theory for the relaxation times over 
viscosities, e.g., and express them in terms of Jmn's. 
Here we note that, as is the case for ordinary hydrody- 
namics, these coefficients depend on space-time coordi- 
nates through their temperature and chemical potential 
dependences in our formalism. Aside from the differences 
accounted for non-trivialities of multi-component sys- 
tems, the equations have almost the same tensor struc- 
ture as that in our formalism with different coefficients. 
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The similarity in tensor structure is worth-mentioning, 
considering that their formahsm uses the technique based 
on renormalization group theory and thus is different 
from our formahsm. It should be mentioned that this 
formahsm has frame dependence, and the forms of the 
constitutive equations differ from the above ones when 
the particle frame is employed. 



3. Constitutive Equations from Methods of I4 Moments 

The second order single-component constitutive equa- 
tions in the original Israel-Stewart formalism Q are 



n = -C9^< - c/3oi?n 



(84) 



Po T 



KTl3iAf"'Dq^ 



eo 



(85) 



(86) 



Note that the metric used in the above paper is opposite 
to ours, which is the source of the negative sign before 77 in 
the first term in Eq. (|86p . Here single conserved current 
systems with binary collisions are considered. Compared 
with our formalism, they discards the second order terms 
with respect to the first order thermodynamic forces, i.e., 
D^, D^, S/f^uf^, V^^, V^f;, and V(^u^.), stating they 
would be small. It is worth-mentioning, however, that 
the terms proportional to acceleration Dm^ are kept in 
the equations. 

Several post Israel-Stewart second-order terms are 
found for single component systems in Ref. [isj by phe- 
nomcnologically expressing the entropy currents in terms 
of the dissipative currents up to the second order. The 
constitutive equations are, in the particle frame, 



n 



n 



(87) 



(88) 



(89) 



The new terms here are the ones proportional to IIV^u^, 
q^Vi/u" and ■n^'^V pu'' . Terms which involve acceleration 
Du^^ do not appear in the second order terms. It is also 
mentioned in the paper that when the kinetic approach 
is employed, vorticity terms TqUj^^'^q,^ and Tt^tt^'^lo"'^^ ap- 
pear in Eq. ([55]) and in Eq. respectively. Note 
that due to the ambiguities of the phenomenological ap- 
proaches mentioned in Sec. IIII Al bulk pressure 11 does 
not appear in the constitutive equations of shear stress 
tensor tt'^'^, and vice versa. 

More second order terms are reported in Ref. [2l| in 
the framework of Israel-Stewart theory with Grad's 14- 
moment method. Actually when our formalism is re- 
duced to single conserved current systems and the non- 
trivialities of the multi-component systems are omitted, 
the third moment equations from the two formalisms 
should become equivalent except for transport coeffi- 
cients because both of them follow Israel-Stewart theory 
consistently. The apparent difference from our formalism 
comes from expansions. They expand first and then 
take the derivative 



- d^{J^''"Pep)~d^{J^''"'f"^ep^), (90) 

whereas we take the derivative first, then expand it as 
shown in Eq. (^5]) . Here ep-^ = ep-y — -^^^^l^ff^-y is a 
traccless tensor. The coefficient e can be identified with 

Tr(Eg-,) 
4 



in single component systems. 

The two expansions should yield the same results up 
to the second order in small quantities, because both of 
them expand the derivative of the distribution as 



daf = - foil ± fo)dayQ - foil ± fo)daSy 

+ /o(l±/o)(l±2/o)Ma2/0- 



(91) 



The first term of the right hand side corresponds to the 
Navier-Stokes limit. The second term is the source of 
the Israel- Stewart second order terms which includes the 
derivatives of the dissipative currents. Note that all the 
terms derived from this second term do not appear in the 
original paper by Israel and Stewart as mentioned ear- 
lier. The third term corresponds to "new" terms that are 
not shown in Ref. These terms cannot be neglected 
because they are also second order terms. 

We again consider the energy frame to compare our 
multi-component results with single component ones, 
i.e., q^' = -7^^^^- According to Ref. 0], the con- 
stitutive equations arc 



H = -CV^it^ - Tni^H 

- C^nV^w^ + mqqi^Du^' - luqdpq^ 



+ A 



(92) 



14 



eo + Po T 



+ A,nnV^^+T,nn-Du^ + ;,nV^n, (93) 



+ 2A^nnV<^u'^\ 



(94) 



which should be identified with our constitutive equa- 
tions by taking into account the ideal hydrodynamic re- 
lations (l80l) and 



\dno 



"0 



eo + Po T' 



(95) 
(96) 



and the identities ([56)) and ((58)) . 

A generalization of Israel-Stewart theory to a relativis- 
tic gas mixture is investigated in Ref . [ll| . Their formal- 
ism is different from ours in several aspects. The systems 
with no particle creation or annihilation are considered 
in that paper, i.e., the numbers of each particle species 
are conserved, while we consider the systems with conser- 
vations based on quantum numbers. Their derivation of 
the equations of motion explicitly depends on the Boltz- 
mann equation, and the dissipative currents are split for 
each components. On the other hand, our formalism aims 
derivation of the dissipative hydrodynamic equations and 
dissipative currents cannot be split. In the single compo- 
nent limit that formalism reduces to the original Israel- 
Stewart formalism without acceleration Du^^ and vortic- 
ity UJ^"'. 



and unknowns would be lost and the system could be de- 
scribed only when sets of assumptions are made, as one 
can see in the case of single component systems with one 
conserved current and no chemical interaction. Secondly, 
in systems with conserved currents, additional moment 
equations, which are the second moments of the distribu- 
tion /* with conserved charges as weight factors, should 
be introduced. This allows us to describe the system 
completely. Thirdly, Grad's moment method for the de- 
termination of the distortion of the distribution, which 
is necessary for the derivation of the second order con- 
stitutive equations, should also be modified to match the 
number of equations and that of unknowns. We gen- 
eralized the moment method so that the resulting con- 
stitutive equations satisfy Onsager reciprocal relations. 
Fourthly, the law of increasing entropy requires the cor- 
rection tensor to the distribution e^^ to have finite trace 
in multi-component systems. 

We explicitly calculated the entropy production in ki- 
netic theory and made clear that the extended Israel- 
Stewart theory with the moment method indeed satis- 
fies in the first order limit the Onsager reciprocal rela- 
tions which demand the transport coefficient matrices to 
be symmetric. It is important to note that all thermo- 
dynamic forces of the same tensor order appear in the 
constitutive equation for a dissipative current. For ex- 
ample, the charge diffusion due to the spatial gradient in 
temperature is known as Soret effect and energy dissipa- 
tion due to the spatial gradient in chemical potential as 
Dufour effect. We further investigated other formalisms, 
and the phenomenological approaches are found to be un- 
suitable for the derivation of the second order equations 
because they include ambiguities in associating second 
order terms to the constitutive equations. Our second 
order multi-component equations mostly agree with the 
equations of other formalisms in the single conserved cur- 
rent limit, except for the transport coefficients and for the 
presence of independent equations of and Vj in our 
formalism. 
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Appendix A: First Order Dissipative 
Hydrodynamics for Multi-Component Systems with 
Multiple Conserved Currents 

We briefiy review first order dissipative hydrodynam- 
ics for relativistic multi-component systems with multiple 
conserved currents with emphasis on the Onsager recip- 
rocal relation. The entropy current is expanded in terms 
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of 5f = f - /o as 



up to the first order. Then the entropy production is 



SeD- - n-V,,u'' + W^'(s/u- + -Du,, 



- E^/^v,^. 



T 



(A2) 



Here we have used the energy-momentum conservation 
and the charge number conservations for the first equal- 
ity. Then one finds that the dissipative currents should 
be expressed as 



1 



T 



ECnfaj-P^, 



+ Ew.v^^, 

<5"-J = - E CstuStikD^ + pKjnT^Vp 



if 



T 



(A3) 
(A4) 

(A5) 
(A6) 

(A7) 

(A8) 



in linear response theory. Q — is bulk viscosity, 
Kw = Kww energy conductivity, rj = shear viscos- 
ity and KVj = KVjVj charge conductivity of the J-th 
conserved current. Note that different thermodynamic 
forces of the same tensor order are allowed in constitu- 
tive equations, such as the gradient of ^ in the con- 
stitutive equation for W^. For such "cross terms", the 



transport coefficients should satisfy Onsager reciprocal 
relations [s^, e.g., kwVj = i^VjW and have to be so 
chosen that semi-positive definiteness of the entropy pro- 
duction is preserved. 

One might think each term in Eq. (jA2p should be ex- 
pressed as a quadratic form in dissipative currents to 
obey the second law of thermodynamics, e.g., oc 
(V^^ -I- j:DUfj,) and some of the cross terms seemingly 
violate it when Onsager reciprocal relation is considered. 
However we can choose the coefficients so that the law of 
increasing entropy is actually satisfied. If we denote the 
dissipative currents of the same tensor order as Jp and 
the corresponding thermodynamic forces as Xp, then the 
entropy production for this tensor order is 



E/ "^p-^p- 



(A9) 



The dissipative currents in linear response theory are ex- 
pressed as Jp = CpqXq where Cpq = Cqp is an ele- 
ment of the transport coefficient matrix. Thus we have 



dpS^ — ^ ^ XpCpqXq 



(AlO) 



p,g 



On the other hand, one can diagonalize any symmetric 
matrix with a certain orthogonal matrix = as 
Cpq = Yl,r,s PpriCr^rs)Psq- Here One requires the diago- 
nal elements must be positive. Then 

5pS^ = XpPp^{C'^6rs)PsqXq 

= J2X'^{C'Jrs)X: 

r.s 

= 5]c;x;2>o, (All) 



is obtained where X'p = Y^qPpqXq. It means that if 
the appropriate linear combinations of the thermody- 
namic forces are taken, the entropy production can be 
expressed in quadratic forms, which obviously satisfies 
its semi-positive definiteness. 

It is important to mention that the constitutive equa- 
tions for 5e and 5nj have no meanings when the stability 
conditions are employed because they require 5e = 5nj ^ 
0. The constitutive equations are in the form of linear 
response to the thermodynamic force such as V^m^ or 
V^^m''^ Note here that first order theory is independent 
of the specific form of the distribution . 

These first order expressions are acausal and unstable 
3 because they allow instantaneous propagation, and 
one has to consider the second order corrections to 
to obtain a causal theory. The expansion of the entropy 
current at the second order yields 



-E 

2 ^ 



(A12) 
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which is negative when contracted with the flow m^, rep- higher order contributions, expressed as follows: 
resenting the fact that system is in maximum entropy 
state. This is the thermodynamic stability condition at 
the second order. Its derivative di^^Ss'^^^ obviously requires 

the information of Sf^. This suggests that, unlike first 
order theory, second order theory in general depends on 
how the distribution is estimated. 



Appendix B: Second Order Distortion of 
Distribution 



The expansion of the phase space distribution /' — 
[exp (y*) T 1]""'^ around the local equilibrium distribution 
/q up to the second order yields the distortion of the 
distribution 



sf = -m ± fM + 3/0(1 ± r,){i ± 2r,)5f'. m 



We estimate in the extended moment method as 
<^y* = PiY^jQi^i +PiPi£fj,vy and the distortion up to 
the second order is expressed as 



6f = 



,/ 

I foil ± /o)(i ± 2/o)(pfpr E if 'if 4^^ 



J,K 



In principle, the 10-|-4iV unknowns are again determined 
by matching the macroscopic dissipativc currents with 
the ones calculated in kinetic theory. Unlike in first or- 
der theory, though, the distortion involves higher order 
contributions in terms of the dissipative currents. The 
resultant and e^y would be, by omitting the third and 



,7 ^ 

Snll + BseSe + E BsnjSnj j Uf^u^ 
/^nnll^ + /SseSeSe'^ + E I^SruSuK SnjSriK 

J.K 

J J 

PwwWPWp + E Pwv,WPVp' 
J 

K ' 
/^nnll^ + fi&e&e&(? + E ~P&ri.,i,nK SnjSriK 

Pii&e^Se + E /3n5n J 11(571 ,/ + E i^SeSnjSeSnj 
J J 

PwwWWp + E PwvjWVp' 
J 

J,K ' 

2/3nwniy(^77,)-f 2E/3 

2^5eiv<5eiy(p7t^) -I- 2 ^Sf,y,8eV^^Uy-) 
J 

2 E ^&njW^^{iLUv) + 2 E ^&^KiVK^njV^^Uy-) 

J J,K 

J 

J 

PwwW^Wy + E PwvjW^Vj 
J 

E PvjV^V^V,^ + ^..tt/tt^p, (B3) 



J,K 
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K 

K K 
K 

K,L ' 
X K 

+ E '^/ey.^e^ + E ^U<vMKVi 

K K,L 
K 

where J's and /3's are the second order prefactors which 
can be calculated in kinetic theory. These new terms 
might add higher order contributions on the left hand 
side of the second order constitutive equations (|5T|) - ([54)) 
which can be absorbed into the right hand side by using 
the Navier-Stokes expressions for the dissipative currents 
n, Wt", and nt^". 

Appendix C: Multi-Component Systems with Single 
Conserved Current 

If the system has only one conserved current, then 
the number of dissipative currents match that of un- 



knowns within the framework of the conventional Grad's 
14-moment method even if inelastic collisions are present. 
Then we have 



Sy' = pre,,+pfp,rep., (CI) 



instead of Eq. (|19p . The expression for the entropy pro- 
duction is modified to 



= e^..^^/^"^" > 0, (C2) 



which means the second moment charge-weighted mo- 
ment equations 9q/j" — Yj cannot be constrained from 
the second law of thermodynamics. If we assume that 
Yj is a macroscopic quantity which can be expressed as 
a linear combination of all the dissipative currents, we 
obtain the same constitutive equations as Eqs. (|5ip - ([54)) 
with different transport coefficients. For the scalar, the 
vector and the tensor dissipative currents, Onsager's re- 
ciprocal relations can be satisfied if we are allowed to 
adjust the transport coefficient matrices C^'s and if the 
semi-positive definiteness of the resulting transport coef- 
ficient matrices Ai is preserved. Since this formalism is 
not equivalent to the multiple conserved current theory 
we developed, it should be investigated which formalism 
we should follow to discuss the single conserved current 
systems. 
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